Nonrelativistic QED approach to the bound-electron g factor 
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Within a systematic approach based on nonrelativistic quantum electrodynamics (NRQED), we derive the 
one-loop self-energy correction of order a(Za)^ to the bound-electron g factor In combination with numerical 
data, this analytic result improves theoretical predictions for the self-energy correction for carbon and oxygen 
by an order of magnitude. Basing on one-loop calculations, we obtain the logarithmic two-loop contribution 
of order a^{Za)'^ ln[{Za)^'^] and the dominant part of the corresponding constant term. The results obtained 
improve the accuracy of the theoretical predictions for the IS' bound-electron g factor and influence the value 
of the electron mass determined from g factor measurements. 
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There has been significant progress in experimental in- 
vestigations of the bound-electron g factor during recent 
years 1 1, 2]. The g factor value in these measurements is ob- 
tained from the ratio of the electronic Larmor precession fre- 
quency o^L and the cyclotron frequency of the ion in the trap 
ujc, according to 
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where e is the elementary charge, q is the charge of the ion, 
TOion is the ion mass, and m is the electron mass. The present 
accuracy of the experimental results for carbon and oxygen is 
already below the level of 1 part per billion and is likely to be 
improved in the near future. In order to match the experimen- 
tal precision achieved for H-Uke ion^ accurate calculations 
of the one-loop self-energy jslMl^la]^ vacuum-polarization 
1 4, 7], and nuclear-recoil [§(, ]^ corrections have been per- 
formed during the last decade. An important result of these 
studies is the possibility to extract the electron mass from the 
experimental value for lul/^c according t o Eg . Q. Such a 
determination was presented in Refs. fffsi fioll based on the 
experimental result for H-like carbon fv\ and oxygen. It pro- 
vided an improvement of the accuracy of the electron mass by 
a factor of 4, as compared to the previous value based on mea- 
surements involving protons and electrons in Penning traps 
rUl . As a result, the presently recommended value for the 
electron mass 11211 is derived mainly from the bound-electron 
g factor. 

The current uncertainty of the theoretical values for the g 
factor in H-like carbon and oxygen (as well as for other low- 
and medium-Z ions) originates predominately from the two- 
loop binding QED correction. This correction is presently 
known 11311 only to its leading order in Za, namely a^{Za)^, 
the corresponding contribution being of pure kinematical ori- 
gin. The next correction enters in order c?{Zct)'^ and results 
from non-trivial two-loop binding QED effects that have not 
been addressed theoretically up to now. 

The goal of the present investigation is to formulate an ap- 
proach for the systematic derivation of higher-order QED cor- 
rections to the g factor of a bound electron. Applicability 



of this approach is demonstrated for the one-loop self-energy 
correction, for which direct numerical calculations to all or- 
ders in Za are available. The analytical contribution derived 
here is used in order to reduce the uncertainty of the numerical 
results for carbon and oxygen by an order of magnitude, this 
being the second largest error of the corresponding theoretical 
values for the g factor [5J. The developed method is then ap- 
plied to the evaluation of the two-loop self-energy correction. 
We derive the complete result for the logarithmic contribution 
of order o? (Z a)^ \i\\{Z a)^'^\ and a large part of the corre- 
sponding constant term. The results obtained improve the ac- 
curacy of the theoretical values for the 15 bound-electron g 
factor in carbon and oxygen and influence the electron-mass 
determination derived from the corresponding experimental 
results. 

The computational approach is based on the Dirac- 
Coulomb Hamiltonian that is modified by the presence of the 
free-electron form factors F\ and HUl^ 

H = a- [p-eFi(A)A] +/3TO + eFi(A)An 
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where A denotes the Laplacian operator. This Hamiltonian 
accounts only partly for the interaction of the electron with 
high-frequency photons. The remaining high-energy contri- 
bution can be represented as a local type of effective interac- 
tion. It is obtained by matching the low-energy limit of scat- 
tering amplitudes derived from the Hamiltonian (|2} and from 
full QED in a way that will be discussed below. 

In the calculation of electron self-energy corrections, it is 
often convenient to evaluate contributions due to the high- and 
the low-energy virtual photons separately. The separation is 
achieved by introducing a certain cut-off parameter, which in 
our calculation is chosen to be the photon mass. The low- 
energy part is then identified as the difference of the contri- 
butions with the massless and the massive photons. In order 
to calculate it, we perform the Foldy-Wouthuysen transfor- 
mation of the Hamiltonian iJ followed by the Power-Zienau 
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transformation, as described in Ref. lll4ll . 

Hpz = ^ er-E- — {l + K)(T-B 
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Here, £; = £; 



B = B\ 



and K 



,,^0' - - -ir=0' - - -^^2(0). In the 
Hamiltonian Hpz we neglected the spin independent terms 
and the A-dependence of the form factors. Moreover, the 
terms with k will be needed only in the two-loop calculation. 

We consider now the one-loop self-energy contribution to 
the bound-electron g factor. It is represented by the sum of 3 
parts, Sg^^^ = g[^'' + (jj^' + g'^^^ . The first part comes from 
the free-electron form factors in the Hamiltonian H, the sec- 
ond part is due to an additional term that matches scattering 
amplitudes, and the third part is a low-energy-photon contri- 
bution that originates from the Hamiltonian Hpz, and is very 
similar to the Bethe logarithm in the hydrogen Lamb shift. All 
these parts are calculated in the following. 

Form-factor contribution. — We evaluate this part by sepa- 
rating the Hamiltonian (|2j into the unperturbed (Dirac) Hamil- 
tonian and the interaction part and applying the standard 
Rayleigh-Schrodinger perturbation theory. Taking into ac- 
count that only the first 2 terms in the A expansion of form 
factors contribute to the order of interest, we write the inter- 
action Hamiltonian as 



SH 



e K 
2m 



ij-E-f3'E-B 



+eFi{0)AAo + —F^{0)A{j-E), (4) 
2 m 

where the slope of the form factors are known to be F[ (0) = 
77 (-1/4 - 2/3 In^) and F^(0) = a/(12 7r), where ^ is the 
ratio of the photon mass to the electron mass. Applying per- 
turbation theory in first and second orders and separating con- 
tributions linear in the magnetic field, we obtain the correction 
to the g factor of an nS state. 



9i 



(1) _ 



1 



7 5 16 , 

- H 1 In u 

6 24n 3 ^ 



6 

' (5) 

Spin-dependent scattering amplitude. — It represents a 
high-energy contribution which goes beyond the on-shell 
form-factor treatment. Here we only sketch the idea of the 
derivation; the details will be presented elsewhere. We first in- 
troduce the skeleton amplitude of the free electron scattering 
off both the Coulomb and the magnetic external field. Then 



we add an electron self-energy loop inserted into the skele- 
ton diagram in all possible ways. Infrared divergences present 
in loop-momentum integration are regularized by employing 
the photon mass. Next, we subtract the skeleton amplitude 
with vertices modified by the electron form factors Fi and 
F2, expand in all external momenta, and keep terms of the 
third power in external momenta only. The resulting ampli- 
tude can be represented by the following effective-interaction 
Hamiltonian (B = const.) 



SH2 = 



a 



— cTiBt djEj + — Irifi Oi Bj di Ej 
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^2") The corresponding correction to the energy reads 



(-Za^) S^r)ea-B[l + -lnf,j\<j)), (7) 



which induces the following contribution to the g factor 
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Low-energy part. — This contribution is induced by the vir- 
tual photon of low energy. We first recall the expression for 
the low-energy part of the Lamb shift. 
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where Hs is the Schrodinger Hamiltonian. In the above, we 
assume the implicit difference between massless and massive 
photons and keep only the terms that do not vanish when 
/i 0. In practice, one calculates this with a cut-off k < e 
and later performs the replacement In 2e In /i + s/e ITsll . 
We mention that this replacement is not unique and its spe- 
cific form depends on the actual integrand. Eq. (|9} coin- 
cides with the standard definition of the Bethe logarithm In ko, 
which has the expHcit values In fco(lS') = 2.984 128 555 and 
lTiko{2S) = 2.811 769 893. 

Here, we are interested in all possible relativistic correc- 
tions to Eq. (|9j induced by the Hamiltonian Hpz which are 
linear in the B field. There are 6 such corrections presented 
in Table I. The terms involving E and djB^ represent correc- 
tions to the vertex (— e r ■ E), and the others yield corrections 
to H, E, and (p. The results listed in the third column of Table 
I involve the standard Bethe logarithm In fco and its modifica- 
tion In fcs defined as 
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We calculate In k^ numerically, using a finite-difference rep- 
resentation of the Schrodinger Hamiltonian Hs, and obtain 
the following results for the IS" and 25 states: lnfc3(lS') = 



3 



3.272 806 545 and ln/c3(2S') = 3.546 018 666. Finally, we 
present the sum of all 6 contributions from Table I as 
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Two-loop contribution. — The complete two-loop self- 
energy of order {ZaY consists of contributions related to the 
electron form factors, to the two-photon scattering amplitude, 
and to the low-energy part. In the present investigation we 
derive an expression for the low-energy part only. This gives 
the complete result for the logarithmic contribution of rela- 
tive order (Za)^. Moreover, we observe that in the one-loop 
case, the low-energy part yields the dominating contribution 
of about 75% of the constant term. Arguably, this is a general 
feature of all radiative corrections, another example being the 
hydrogen Lamb shift. We thus assume that also in the two- 
loop case, the low-energy part provides the dominant contri- 
bution to the constant term. A derivation of the remaining 
two-loop contributions can in principle be carried out along 
the lines presented above. 

Let us now identify the two-loop low-energy correction. 
When both photons are of low energy, the B-dependent part 
is of a higher order and thus negligible. The contribution of 
interest comes when only one photon is of low energy. The 
second photon effectively modifies the vertex, and only the 
part with the anomalous magnetic moment is relevant, as the 
slope of the form factors contributes to higher orders. There 
are two equivalent contributions obtained by interchanging the 
photons, which results in an additional factor of 2. The effec- 
tive Hamiltonian that accounts for the anomalous magnetic 
moment is given by Eq. (|3}- The calculation of the Bethe- 
logarithmic contributions is the same as for the one-loop case, 
but involves different overall factors for each term, which are 
listed in the fourth column of Table I. Using the one-loop re- 
sults, we obtain for the sum of all low-energy contributions. 
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where the numerical value of the constant term is a^Q^ (15) = 
-18.477948 664(1) and 40^(25) = -19.781820 939(1). 
The term with In p. is canceled by corresponding contributions 
coming from the slope of the form factors and the two-loop 
scattering amplitude. 

Results and discussion. — We first summarize our calcula- 
tion for the one-loop self-energy correction. The total analytic 
result for an nS state is 
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where the remainder function Gn{Z) incorporates all contri- 
butions of higher orders in Za, and the numerical value of the 



constant term in order {Za^ is a^g^ = -10.236 524 318(1) 
for the IS state and a^J,^ = -10.707 715 607(1) for the 25" 
state. The first two terms in Eq. ( I13> are well known; the first 
one is the famous Schwinger correction and the second was 
derived previously for the IS state in Ref. fl3ll . 

By subtracting all known terms of the Za expansion in 
Eq. ( I13> from numerical data one can isolate the one- 

loop self-energy remainder Gn{Z) and improve its numerical 
accuracy for carbon and oxygen by extrapolating results for 
higher values of Z. The higher-order contribution extracted 
directly from numerical results of Ref. |5] reads 6*1(6) = 
22.19(24) andGi(8) = 21.86(6). An extrapolation of numer- 
ical data |5] for Z > 8 yields the results for the self-energy 
remainder Gi(6) = 22.160(10) and Gi(8) = 21.859(4), 
which are significantly more accurate. In addition, we obtain 
the following result for the total contribution of order [Za)^: 
Gi(0) = 23.0. 

The result for the two-loop self-energy contribution is given 
by Eq. Mil . We estimate the uncertainty due to uncalculated 

(2) f2) 

parts g\ ' and g\ ' as 30% of the constant term. Explicitly, the 
two-loop self-energy correction for the IS* state is 

<5g(2) = {^\zaf ln[(Za)-2] - 18.5(5.5)} . (14) 

We now turn to the experimental consequences of our cal- 
culation. A previous compilation of theoretical contributions 
to the 15 bound-electron g factor was given in Ref. d. In 
the present work we modify it in several ways, with the corre- 
sponding contributions listed in Table HII {i) we employ the 
new, more accurate results for the one-loop self-energy re- 
mainder; { ii) we use the analytic result of Ref. f?] for the lead- 
ing term of the Za expansion of the first-order magnetic-loop 
vacuum-polarization correction; (Hi) we include the leading 
part of the two-loop self-energy correction of order a'^{Za)^ 
obtained in this work [Eq. (I14H . We assume that the uncer- 
tainty due to other uncalculated two-loop corrections is ab- 
sorbed into the error bars of the constant term in Eq. J14t . 

As compared to the previous compilation [5], the accuracy 
of the theoretical value for carbon is improved by a factor of 
3. In case of oxygen, only a small improvement of accuracy 
is achieved, but the theoretical value is shifted slightly outside 
of the error bars given in Ref. 0]. The described modification 
of the theoretical predictions for the bound-electron g factor 
influences the electron-mass values derived from the experi- 
ments on carbon QJ ^nd oxygen Following Refs. o fioll 
and using the g factor values from Tablellll we obtain the fol- 
lowing results for the electron mass, 

mC^&+) = 0.000 548 579 909 41 (29)(3)u, (15) 
mC^O'^+) = 0.000 548 579 909 87 (41)(10)u, (16) 

where the first uncertainty originates from the experimental 
value for the ratio ujl/uJc, and the second error comes from 
the theoretical result for the bound-electron g factor. 

In summary, we have presented an approach for a system- 
atic derivation of higher-order QED corrections to the g fac- 
tor of a bound electron. We obtained the complete result for 
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TABLE I: Breakdown of the low-energy contribution to the bound-electron g factor. 
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TABLE IL Individual contributions to the Is bound-electron g factor, 1/a from IT3l is 137.035 999 11(46). 

12q5+ 16q7+ 

Dirac value (point) 1.998 721354 39(1) 1.997 726 003 06 (2) 

Finite nuclear size 0.000 000 000 41 0.000 000 00155 

FreeQED, ~ (q/tt) 0.002 322 819 47(1) 0.002 322 819 47(1) 

Binding SE, ~ (q/tt) 0.000 000 852 97 0.000 001 622 67 (1) 

Binding VP, ~ (q/tt) -0.000 000 008 51 -0.000 000 026 37 (1) 

FreeQED, ~ (Q/7r)2--(a/7r)* -0.000 003 51510 -0.000 003 51510 

Binding QED, ~ (Q/7r)2(ZQ)2 -0.000 000 001 13 -0.000 000 002 01 

BindingQED, ~ (a/7rf (Za)'' 0.000 000 000 41 (11) 0.000 000 00106(35) 

Recoil 0.000 000 08763 0.000 000116 97 

Total 2.001041590 52(11) 2.000 047 02128(35) 



the one-loop self-energy correction of order a{Za)'^. The 
derived contribution is in excellent agreement with the pre- 
vious numerical calculation. The developed approach was 
then applied to the most problematic two-loop self-energy cor- 
rection. We obtained the logarithmic contribution to order 
a'^{Za)'^ \n{Za)^^ and the dominant part of the correspond- 
ing constant term. As a result, we improved the accuracy of 
the theoretical predictions for the IS* bound-electron g factor 
for carbon and oxygen and presented new values for the elec- 
tron mass derived from the corresponding measurements. 
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